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Abstract. Despite a wealth of experimental data for high-P^ processes in heavy- ion 
collisions, discriminating between different models of hard parton-medium interactions 
has been difficult. A key reason is that the pQCD parton spectrum at RHIC is 
falling so steeply that distinguishing even a moderate shift in parton energy from 
complete parton absorption is essentially impossible. In essence, energy loss models 
are effectively only probed in the vicinity of zero energy loss and, as a result, at 
RHIC energies only the pathlength dependence of energy loss offers some discriminating 
power. At LHC however, this is no longer the case: Due to the much flatter shape of 
the parton pr spectra originating from 2.76 A GeV collisions, the available data probe 
much deeper into the model dynamics. A simultaneous fit of the nuclear suppression 
at both RHIC and LHC energies thus has great potential for discriminating between 
various models that yield equally good descriptions of RHIC data alone. 



PACS numbers: 25.75.-q,25.75.Gz 
1. Introduction 

High-JV processes taking place in the background of the medium produced in 
ultrarelativistic heavy-ion (A-A) collisions are a cornerstone of the experimental A- 
A program at the LHC. The aim is to use these processes to do "jet tomography", i.e. 
to study both the short-distance physics of the bulk medium (i.e. its relevant degrees 
of freedom) and the geometry of its expansion. However, so far attempts to extract 
even solid qualitative statements about the nature of parton-medium interaction have 
not been successful. Two main reasons can be identified: 1) For steeply falling primary 
parton spectra, medium-induced shifts in parton energy cannot be distinguished from 
parton absorption, hence observables lose sensitivity to model details; 2) Computed 
observable quantities depend both on assumptions made about the bulk medium 
evolution and the parton-medium interaction. In this work, we propose to resolve 
the second ambiguity by a systematic investigation of multiple observables, while we 
demonstrate that the first is significantly lessened for the harder LHC parton kinematics. 
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2. Dependence on medium modelling 

We test different combinations of medium evolution and parton-medium interaction 
models against a large body of high Pj- observables. In particular, for the medium 
evolution we use a 3+ld ideal p], a 2+ld ideal [2] [3] and a 2+ld viscous hydro code 
[I] with both CGC and Glauber initial conditions. On the parton-medium interaction 
side, we test a radiative energy loss model [5], a parametrized [6] and a Monte-Carlo 
(MC) model [7] for incoherent energy loss, a strong-coupling phenomenological model 
based on AdS / CFT ideas [8] and the MC in- medium shower code YaJEM jH [10] with 
its variant YaJEM-D pj] which introduces an explicit pathlength/energy dependence 
into the minimum virtuality scale down to which the shower is evolved. 
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Figure 1. Left panel: Raa{Pt) for 30-40% cental 200 AGeV Au-Au collisions for in 
plane (solid) and out of plane (dashed) emission computed for the same energy loss 
model (ASW) with different hydrodynamical backgrounds, compared with PHENIX 
data [12]. Right: Raa{4>) at Pt — 10 GeV for smooth, initial-state averaged and 
event- by- event fluctuating hydrodynamics for different fluctuation size scale. 



In a first run, we pose the question to what degree the underlying medium model 
is able to influence high Pt observables. In Fig. [T] we present an example of results 
from a systematic investigation of both the influence of smooth hydrodynamical models 
[T3] and event-by-event hydrodynamics [3] with initial state density fluctuations [14]. 
Summarizing the result, we find that the medium model has a considerable (factor ~ 2) 
influence on observables such as Vi at high Pt or extracted parameters such as the 
transport coefficient q. Chiefly responsible is the location of the freeze-out hypersurface 
- the agreement with data in general improves if the freeze-out hypersurface is large, 
but noticeable effects are also caused by the initialization time or the presence/absence 
of viscosity. Fluctuations in the hydrodynamical initial state play a minor role for 
extracted parameters (~ 20%), due to a cancellation of competing effects (see [14"]). 
but for non-central collisions the cancellation is incomplete, leading to a decrease in 
suppression for small fluctuation size scale. If hard probes are used to constrain the 
fluctuation size, a scale of ~ 0.8 fm is preferred. 
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3. Pathlength dependence 

In a second run, we test the pathlength dependence of the available parton-medium 
interaction models against the data [HI [13]. In a static medium with constant density, 
we expect incoherent processes to scale with pathlength L (elastic), radiative energy loss 
(ASW) with L 2 due to coherence effects and the strong coupling model (AdS) with L 3 
due to the drag-force like interaction of the virtual gluon cloud with the medium. The 
shower code YaJEM is known to have in principle an L 2 dependence due to coherence, 
but which effectively reduces to L by finite-energy corrections, whereas YaJEM-D has a 
complicated non-linear pathlength dependence. In an evolving hydrodynamical medium, 
the pathlength dependence is effectively much more complicated due to effects like 
spatial inhomogeneities in the medium, longitudinal and transverse flow and viscous 
entropy production. We find that the data allow to unambiguously rule out linear 
pathlength dependence, leading to the conclusions that a large component (> 10%) 
of elastic energy loss is not favoured by the data and that finite-energy corrections to 
coherence arguments need to be taken seriously. The other models we tested (ASW, AdS 
and YaJEM-D) remain viable with the data, although in each case only in combination 
with a particular hydrodynamical evolution model. 

4. Extrapolation to large y/s 

Using the EKRT saturation model, we can extrapolate one of our default hydro runs to 
LHC energies [15] and thus significantly reduce the uncertainty in the hydrodynamical 
modelling, as the well-defined extrapolation procedure allows to compare results for 'the 
same' hydrodynamics at different y/s. 
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Figure 2. Left: The nuclear suppression factor Raa extrapolated from best fits at 
RHIC to central Pb-Pb collisions at LHC for various models (see text) compared with 
ALICE data [16]. Right: Raa for non-central collisions in-plane and out of plane based 
on a best fit to central LHC data. 

Fig. [2] shows results for Raa m central and non-central collisions. As expected, the 
sensitivity to model details in the Pt dependence of the nuclear suppression factor is 
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found to be much larger at LHC than at RHIC. Combined with pathlength-dependent 
observables such as Raa{4>) or the dihadron suppression factor Iaa, precision LHC data 
thus have a high potential for distinguishing between different models. 

We can quantify the quality of extrapolation in yfs by the factor R which quantifies 
the differences of parton-medium interaction parameters for a best fit to RHIC and LHC 
data, using the same hydrodynamical models. R = 1 indicates an extrapolation without 
any tuning. We find that the shower codes YaJEM-D (R = 0.92) and YaJEM (R = 0.61) 
extrapolate reasonably well, whereas the radiative energy loss scenario ASW {R = 0.47) 
is not favoured by the data, and a strongly coupled scenario AdS is strongly disfavoured 
with R = 0.31. The latter effect can be readily understood by a dimensional analysis - 
if the pathlength dependence is L 3 , then for dimensional reasons the model must probe 
the medium temperature as T 4 , i.e. the model responds much more strongly than all 
others to the higher initial medium density at the LHC, which leads to overquenching. 
Given this finding, there is currently no reason to assume that the data would prefer a 
strongly coupled over a perturbative scenario of parton-medium interaction. 

Combining the constraints from pathlength dependence and \fs extrapolation, 
assuming no systematic uncertainty on the published LHC data [16] , out of the models 
tested here only YaJEM-D together with a hydrodynamics similar to the 3+ Id code 
remains a viable description. We take this as a strong indication that systematic studies 
along the lines discussed here are indeed a suitable tool to do jet tomography. 
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